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Abstract 

Star formation rates on the galactic scale are described phenomenologically by two 
distinct relationships, as emphasized recently by Elmegreen (2002). The first of 
these is the Schmidt law, which is a power-law relation between the star formation 
rate SFR and the column density S. The other relationship is that there is a cutoff 
in the gas density below which star formation shuts off. 

The purpose of this paper is to argue that 1) these two relationships can be 
accommodated by a single functional form of the Schmidt law, and 2) this functional 
form is motivated by the hypothesis that star formation is a critical phenomenon, 
and that as a corollary, 3) the existence of a sharp cutoff may thus be an emergent 
property of galaxies, as was argued by Seiden (1983), as opposed to the classical 
view that this cutoff is due to an instability criterion. 
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1 Introduction 



The Schmidt law is a power-law relationship between the star formation rate 
SFR and the gas column density E, introduced by Schmidt (1959), namely 

SFRocS™ (1) 



for some exponent n. Schmidt estimated n to be close to 2; modern mea- 
surements of the Schmidt index tend to give lower values for n, such as 1.4 
(Kennicutt, 1998; Heyer et a!., 2004), although the index depends upon what 
material is included (neutral gas, ionized gas, stars) in S. Note that Wong 
& Blitz (2002) indicate that the most reliable tracer is the surface density of 
molecular hydrogen; even so, there is still significant variance in inferred values 
of n, in part due to observational and modeling uncertainties, but presumably 
in part due as well to real variation in Nature. 
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The other relationship is that t here is a cutoff in the gas density below which 
star formation shuts off. Again, E lme green (2002) emphasizes that this cutoff 



is either interpreted phenomenologically as a fixed gas density, or a density 
determined by a stability criterio n (e.g. the Toomre Q criterion). The history 



of this cutoff is explored nicely in iHunter et all ([199 



The basic Schmidt Law of eq. (1) has many variants, which we will not de- 
scribe here. The consistent notion is that, averaged over some suitably large 
lengthscale, the star formation rate in a galaxy increases with the gas density, 
and that it does so as a power law. What is surprising is how well this pre- 
scription seems to work, over a wide range of densities, from spiral galaxies 
up through starburst galaxies. Given the simplicity of eq.(l), it is perhaps 
not too surprising that a wide variety of physical mechanisms have been sug- 
gested as the root of this behavior. Furthermore, various modifications to this 
functional form have been suggested. The original prescription, as well as sub- 
sequent modifications, all share the property that the power law is centered 
on the origin. In other words, the two quantities of concern are connected by 
a relationship of the form 

(x-x )=A(y-y f (2) 



with the particular choice x = and y = 0. 

The problem with this is that it is at odds with the second observational fact 
described above, namely the existence of a cutoff. The result is that, as seen in 
observations, star formation falls off of the power-law curve as one approaches 
the cutoff from above. This drop-off of the power-law is also seen in simulations 
that show a cutoff (Kravtsov 2003). 

The existence of a sharp cutoff in star formation below a critical surface density 
in spiral galaxies is shown in the data of Kennicutt (1989) and of Martin 
& Kennicutt (2001), among others (although note that Pohlen et al. (2002) 
argue against a sharp cutoff). Dynamical instability (Toomre, 1964; Quirk, 
1972; Wang & Silk, 1994; Pandey & van de Bruck, 1999) of some sort may 
explain this cutoff, as argued by Kennicutt (1989) and others. On the other 
hand, the results of Wong & Blitz (2002) and Boissier et al. (2004) do not 
offer as strong support for the s tability criter i on int erpretation of the cutoff 
as might be hoped. In addition, IHunter et al | (|l998t) point out that it is not 



clear that the instabilities used to justify theoretically the existence of a cutoff 
critical density E should actually apply directly to the star-formation process. 
Furthermore, recent star formation may actually have taken place in spirals 
beyond the nominal cutoff radius, so that the observed cutoff may be to a 
degree an issue of observational sensitivity (Ferguson et al. 1998). 
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It is interesting to note that star formation behaving as a power-law with 
lower cutoff is also seen in certain extremely simple models of galactic-scale 
star formation, namely self-propagating star formation (SPSF) and stochastic 
self-propagating star formation (SSPSF). Such models will also show a small 
but significant degree of star formation below the cutoff, depending on the 
parameters of the simulation. On the other hand, it is not at all clear that 
star-formation on galactic scales proceeds by self-propagation. For example, 
one current view is that supernovae (SN) are effective pumps of supersonic 
turbulence in potential star-forming clouds, and that such turbulence provides 
the dominant means of support against Jeans collapse and star formation 
(Mac Low & Klessen, 2004), implying that star formation in one locale is 
hindered rather than helped by recent star formation in nearby locales. 

Here we do not advocate SPSF or SSPSF one way or another. However, we 
do take some lessons from the perspectives these models provide, such as the 
point of view that the star formation cutoff may an emergent property of the 
system, and might not be reflective of a classical dynamical instability. Instead, 
as has been argued previously (see below), star formation and cutoff may be 
similar to critical phenomena more familiar to condensed-matter physicists. 
In such critical phenomena, the cutoff is incorporated into the power law, as 
shown below. 



1 . 1 SPSF and SSPSF models 



As motivation to the points of view alluded to above and discussed further 
below, we now outline the basic considerations of stochastic self-propagating 
star formation models. It should well be borne in mind, however, that the phe- 
nomenon of a critical point in global star formation rates does not necessarily 
imply that star formation is indeed a self-propagating phenomenon; rather, 
this specific instance of critical behavior in global star formation models is 
simply offered as motivation for our modified functional form. 

The self-propagating star -formation model (SPSF) was originally proposed by 
Mueller fc Arnett ( 1976| ) as a simple model that offered an explanation of spi- 



ral patterns in disk galaxies as the result of two processes, namely the presence 
of large-scale shear in disk galaxies, and the hypothesis that star formation in 
one lo cation induces star formation in neighboring regions. Gerola &: Seidenl 



(1978) introduced the stochastic self-propagating star-formation model (SSPSF 
model), which modified the model of Mueller & Arnett so that the induce- 
ment of star formation is probabilistic in nature. The SPSF/SSPSF inter- 
pretation of spiral arms may be particularly apropos in the case of so-called 
"flocculent" galaxies, which appear rather spatially disordered, as opposed 
to "grand-design" spirals in which the spiral arms are globally coherent and 
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appear to be dynamical in nature. Both deterministic and stochastic SPSF 
models have been explored more or less continuously in the lite rature since 
the publication of these two papers. Schulman fc Seiden (1986) provide an 



early review of the basics of SPSF models. 



One o f the aspects of the SSPSF model that was noted early on bv lSeiden et al 



( 1979b is that the equilibrium star formation rate as a function of the control 



parameter P st (which parametrizes the probability that star formation is in- 
duced in a given cell given that there is star formation in a neighboring cell) 
exhibits a phase transition: The onset of star formation in SSPSF models as a 
function of P st is a critical phenomenon, in which the star formation rate SFR 
is an order parameter (one of several) and the probability P st is the control 

(C) 

parameter. Below some critical probability P st the star formation rate SFR is 
very low, becoming essentially zero far below the critical point. As one nears 

(C) 

the critical point P st from below, the system undergoes larger and larger 
fluctuations; these fluctuations may also be seen above the critical point in 
the case that the computational domain is very small. For example, the star 
formation may exhibit large bursts between protracted latency periods with 
no star form a tion, as discussed in the theory of d warf galaxies discu ssed in 
Gerola et al (1980) and subsequent papers and by Hirashita ( 2000j ). These 



bursts are initiated by random star-formation seeding events that are con- 
trolled by a second parameter, the spontaneous seeding rate, that is of little 
importance above the critical point so long as it is sufficiently small. 

Far above the critical point, star formation is more or less constant in time. 
Approaching the critical point from above, star formation as a function of 
time becomes more and more erratic, and in the absence of seeding events 
it shuts off entirely below the critical point. Most importantly, the average 
(as a function of ti me and space) star formation rate < SFR >, as noted by 
Seiden et al behaves as 



gFR = / A(P 8t — Pg^)^ if P s t > Pst^; (3) 
[ otherwise, 



which is characteristic behavior for a critical phenomenon (see fig. la). Other 
properties of the system (such as the mean size of star-formation clusters, the 
variance in the star formation rate, etc) also depend upon P st with the same 
functional form, albeit with different critical indices. 

These relationships hold best near the critical point, but computational re- 
strictions limit one from reproducing this behavior arbitrarily close to that 
point: Any real computation is necessarily limited to a finite number of cells, 
and this affects the dynamics of the system very close to the critical point, 
because correlation lengths diverge there. Furthermore, the small but nonzero 
rate of spontaneous star-formation that is an ingredient of most SSPSF models 
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blunts the cusp in the phase transition at the critical point, just as, for ex- 
ample, the random-field Ising model departs from canonical critical behavior 
close to the critical point in that system. Also, behavior of the system de- 
parts from relationship (3) far above the critical point; eventually the system 
saturates and SFR bec omes less dependent upo n P 5t . However, let us keep i n 
mind that, as noted by ISeiden et all (|l979h and ISchulman fc Seiden 



the spiral patterns that were the original motivation for the SSPSF model are 
best produced with P st just above the critical value. 

As a demonstration of the functional form (3), in figure (A.l) we show SFR 
as a function of P st in a very simple SSPSF code we have written. This code 
implements SSPSF in a fixed 2-D (400 x 400) Cartesian lattice with periodic 
boundary conditions. Simulations are followed for 10000 time steps, and the 
initial 1000 time steps are discarded from the analysis to remove transient 
behavior. Cells are in one of four states: ready for star formation, forming 
stars, supernova, and dormant. On each timestep, for a cell that is "ready," 
the number of nearby (above, below, left and right) supernova are added up, 
and multiplied by the parameter P st . If a random number drawn from the 
uniform distribution on [0,1] is less than this number, then star formation is 
initiated on the next step. The cell then goes "supernova" on the subsequent 
step, then it is "dormant," and finally it returns to being "ready." There is 
random spontaneous seeding of star formation inserted at a very low rate 
(roughly one per 20 time steps for the entire mesh in this case). No attempt is 
made to account for the gas consumed by star formation. This is an extremely 
simplified version of SSPSF, used simply to illus trate the form in eq. (3) for 
this system, as shown also in lSeiden et al Here we find the exponent 



(3 to be 0.52, close to the mean-field value of 1/2, when we fit over the range 
0.0 < P < 0.50. At higher values of P the star formation rate begins to 
saturate. 

Just as this cutoff in SSPSF models is called a phase transition, so too is the 
cutoff in observations of star formation referred to as a phase transition. 



1.2 Hypothesis 



Our hypothesis, based on the above notions, then, is that star formation should 
be written as a power-law function not of the gas density, but of the difference 
between the gas density and the critical (or threshold) gas density, in analogy 
to the behavior of critical phenomena: 

SFR oc { ^ ~ ^ n if S - (4) 
1 otherwise, 
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Here, the density S plays the role of the control parameter. We do not address 
here the extensive literature on the various universality classes in the theory of 
critical phenomena, and we make no attempt to connect the exponent (3 from 
SSPSF models to the Schmidt index n. Rather, the point is merely that, within 
the theory of critical phenomena, power laws are always expressed about the 
critical point, and so insofar as star formation may be thought of as a type of 
phase transition, it seems natural to incorporate the cutoff into the Schmidt 
law in the form of eq. (4). 

Note that both eq. (1) with cutoff and eq. (4) have the same number of free 
parameters, namely three: An amplitude (the proportionality constant, not 
written), a power- law index, and a cutoff. 



2 Data Analysis 



To test the relative merit of eq. (1) with cutoff and eq. (4), we fit two sets of 
data to these two functions. The x 2 f° r these fits is found by assuming that 
the uncertainties in the logarithm of the physical quantity on the abscissa and 
the ordinate are equal and constant for all data points. This greatly simplifies 
the fitting process described below. Rather than tabulate the \ 2 f°r the fits, 
we give the ratio of the x 2 f° r the best fits for the two functional forms, since 
this ratio does not depend upon the actual value for the uncertainties of the 
measurements. 



To fit the ordinary Schmidt law, we use established methods for fitting to a 
straight line with uncertainty in both abscissa and ordinate, with the following 
modifications: A line is drawn perpendicular (in log-log space) to the fitted 
line at the cutoff point. Points to the left of this line contribute to the \ 2 of the 
fit by an amount proportional to the distance of these points fro m the cutoff 



point. We then minimize \ 2 using a downhill simplex method ((Press et al 
(1992)). 



To fit the modified Schmidt law, we estimate a relative x 2 by finding the 
closest straight-line distance on the log-log plot between each point and the 
fitted curve. This requires a function minimization for each data point. Again, 
we minimize \ 2 using a downhill simplex method. 

The first set is the data compiled by Hunter et al. (1998), which includes Ha 
surface brightness as a proxy for star formation rate and total (atomic and 
molecular) gas density inferred by linear proportionality to measured HI for 
a range of irregular galaxies. We throw out two galaxies from this data set. 
One, DDO 105, does not extend in radius out past the star formation cutoff. 
The other, DDO 155, has few data points in the star-forming region. 
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Table 1 

Relative Merit of Fits for Data of Hunter et al. (1998); see text for definitions of a, 
n, and c. 







standard 






modified 






object 


a 


n 


c 


a 


n 


c 


^mod/ Xstd 


DDO 154 


26.5926 


5.4241 


0.7776 


30.9764 


1.5926 


0.7546 


1.6148 


DDO 168 


29.8522 


1.6172 


0.5434 


31.201 


0.6682 


0.6965 


0.184 


DDO 50 


26.1368 


5.5392 


0.8159 


30.2791 


2.4847 


0.7864 


1.0009 


IC 1613 


31.009 


7.8061 


0.4186 


17.567 


24.0863 


0.4779 


0.927 


Sex A 


22.1768 


10.8149 


0.5042 


26.4300 


7.5003 


0.4004 


1.0246 



These are dwarf irregular galaxies, and star formation in such galaxies differs 
from star formation in spirals in many respects. Indeed, dwarf galaxies do not 
show a strong correlation between Ha and HI (Brosch et al. 1998). Following 
Wong & Blitz (2002), it may be that much tighter correlations would be 
apparent here if data were available for molecular hydrogen (H2) instead of 
atomic hydrogen. 

The data as plotted here consists of Ha surface brightness in erg/sec/pc 2 
as a measure of star formation rate versus gas column density quantified in 
M /pc 2 . Fitting parameters a,n, and c correspond to the formulas 

log 10 (S(Ha)) = { a + n lo Sio( S g -) if E gas > Eg, ; (5) 
[ — 00 otherwise, 

where c = log 10 (E gas ), in the case of a fit to a standard-form Schmidt law, and 

log 10 (E(Ha)) = { a + n lo Sio( S g as " E ffi) if > Eg) ; (6) 
m v n \ -00 otherwise, v 1 



in the case of a fit to the modified Schmidt law of the form of eq. (4). In 
particular, the quantity n here is thus the usual Schmidt exponent. 

Table 1 shows that the results are equivocal here. For some galaxies, a simple 
power-law fit with cutoff is a better fit, whereas for others our modified power- 
law form appears to be a better fit. Clearly, for some galaxies, both equations 
offer quite poor fits. 



For our second analysis, we choose to reconsider the data of Kennicutt (1989). 
In particular, we analyze the data taken from his figure 8 on page 694 of 
that paper. These data are measured Ha surface brightne ss versus pu b lished 



hydrogen (HI + H 2 ) gas densities, and are not available in Kennicutt ( 1989h 
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Table 2 

Relative Merit of Fits for Data of Kennicutt (1989) 

standard modified 



object 


a 




n 


c 


a 




n 




c 




all 


29.59 


2 


352 


0.3255 


31.36 


1 


085 





5301 


0.606 


M101 


29.28 


2 


944 


0.2533 


30.45 


1 


952 





2422 


0.962 


M51 


29.20 


2 


297 


0.3150 


30.25 


1 


639 





3797 


0.771 


NGC 4254 


30.84 


1 


372 


0.6207 


31.10 


1 


224 





2039 


1.061 


NGC 4303 


31.14 


1 


277 


0.4325 


32.11 





670 





6762 


0.647 


NGC 4321 


30.92 


1 


336 


-0.3723 


31.68 





791 





2028 


0.178 


NGC 4535 


30.03 


2 


124 


0.2068 


31.38 


1 


072 





4679 


0.476 


NGC 6946 


28.67 


3 


027 


0.6132 


31.97 


1 


435 





6533 


0.287 



in tabular form. Therefore, we simply recreate the original data by measuring 
the posi tions of each of th e data points on the figure. A similar plot is also 
found in K ennicutt however this plot does not clearly have individual 



data points marked, and we were unsuccessful at attempting to recreate the 
data plotted there. 

Plots of these fits are included in the appendix. 

It can be seen from these fits that the Schmidt exponent is unusually high in 
several cases. This can be attributed in part to fitting unabashedly the entire 
curve, including downturn that is the subject of this paper. Note that the 
exponents for the modified Schmidt law are systematically lower, as would be 
expected. 

Note as well that the results here are not nearly so equivocal as before; in 
particular, the relative merit of fit to the modified SFR law is markedly better 
here than for the standard form, as demonstrated by the \ 2 ratios for the two 
fits. 



3 Conclusions 



We have offered here an interpretation of the star formation cutoff as an 
emergent prop erty of galactic systems, in the vein of SSPSF models, following 
on the work of Seiden ( 1983| ). This motivates us to consider a modification to 
the Schmidt law, such that star formation is written as a power law around the 
critical density, in analogy to the functional form of quantities near a critical 
point. This new functional form captures the downturn in star formation rates 
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at gas densities near the cutoff point. We find that this functional form is a 
better fit to the data for spiral galaxies, but for irregular galaxies the results 
are ambiguous. 

While this arg uably provides some evidence in support of the hypothesis of 



Seiden I (|1983f ). we note that there are several other interpretations of the 



downturn of star formation rate near the cutoff. 

In particular, some other interpretations include: 1) Star formation occurs 
through a global instability with a critical density. Given a quadratic form 
for the dispersion relation, the growth rate of the instability, as a function 



of density, is proportional to y E gas — Eo- If star formation is proportional 
to the density times the growth rate of the instability, we obtain SFR oc 



E gas y E gas — S , which asymptotes to an exponent of 1.5 and, moreover, ap- 
proaches zero nicely from above. (Note that this is one of a class of relations of 
the form SFR oc E^ as (E gas — E ) n_/3 .) 2) The star formation rate is intimately 
connected with the high-density tail of the density PDF of supersonic turbu- 
lence in the ISM, as suggested by Kravtsov (2003). The functional dependence 
of star formation on density is reflective of this tail, which has a power-law 
form with an upper cutoff, resulting in a power-law dependence of star forma- 
tion on gas density with a lower cutoff. In fact, we show in the appendix that 
our modified functional form appears (by eye) to be a much better fit to the 
simulations of Kravtsov (2003) than the standard Schmidt law. This cannot 
be due entirely to feedback processes (as in SSPSF models) because Kravtsov 
(2003) finds a downturn and cutoff even when feedback is turned off. 3) Star 
formation really does have a power-law dependence with instantaneous cutoff, 
and the appearance of a steepening in the Schmidt index near the cutoff is 
simply due to observational mixing of sub-critical and super-critical regions. 

Thus, our functional form does not discriminate between the cutoff as an 
emergent property of the system, and other hypothesis regarding the origin of 
the cutoff. However, we suggest that the critical phenomenon interpretation 
predicts that a variety of different quantities depend upon the parameter P st 
in the form of eq. (3), and comparison of the behavior of these quantities with 
respect to parameters such as the gas density E may offer further evidence in 
support of the critical-point hypothesis for the origin of a cutoff. The devia- 
tions from a sharp cutoff in star formation that have been observed in extreme 
outer regions of disk galaxies as noted above (Ferguson et al. 1998) actually 
may be evidence for, rather than against, a critical phenomenon interpreta- 
tion, as sytems exhibiting critical phenomena deviate from the form eq. (3) 
both if the system is subject to noise and if the system is small enough not 
to approach the thermodynamic limit. Indeed, as the number of star forma- 
tion sites in galaxies is nowhere near large enough to approach the nominal 
thermodynamic limit, unlike laboratory systems (~ 10 23 particles), it should 
not be surprising that the cutoff is not so sharp. However, accomodation of 
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this fact clearly demands modification of eq. (3) through the introduction of 
another parameter to account for such finite-size effects, and is beyond the 
scope of the work presented here. 

Much of this work was performed while the author was a graduate student at 
the University of Texas at Austin. The author thus wishes to thank the theory 
group in Austin for their support. He also wishes to thank D. Hunter and 
A. Kravtsov for kindly providing results of their observations and simulations, 
respectively. 
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A Fits to Data 



Here we show fits to data and simulations. First, we show results from our 
simple SSPSF implementation (fig. la). Next, we show fits to the data com- 
piled by Hunter et al (1998) for irregular galaxies (fig. 2, c-g) We show fits to 
Ha surface brightness for concentric annuli, used as a proxy for star formation 
rate. The abscissa is total hydrogen (HI + H 2 ) gas density, inferred from mea- 
sured HI emission. We then show fits to the data collected in Kennicutt (1989) 
for Sc spiral galaxies (fig. 2 and 3, h-n). Here the abscissa is the total hydrogen 
(HI + H2) gas density, where atomic and molecular densities were taken from 
previous studies in the literature. Finally, we show fits to the simulated data 
of Kravtsov (2003) (fig. lb). Note that in this final case, the relative x 2 for the 
goodness of fit, i.e. Xmod/Xstd; i s j us t barely less than one, namely it is 0.956. 
However, this measure fails to capture the better "Chi-by-eye" feature of the 
points that they are curved, and that this curve is captured by the modified 
functional form. 
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Fig. A.l. (a) Results and fit of eq. (3) for a simple SSPSF implementation described 
in the text, (b) Fits of the standard Schmidt law, eq. (5), and the modified form, 
eq. (6) to numerical simulation results of Kravtsov (2003). 
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Fig. A. 2. Fits of standard and modified forms (eqns. (5) and (6)) of Schmidt law to 
data compiled by Hunter et al (1998) (c-g) and Kennicutt (1989) (h-n). 
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Fig. A. 3. Fits of standard and modified forms of Schmidt law, continued. 
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This figure "figAl.png" is available in "png" format from: 
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This figure "figA2.png" is available in "png" format from: 
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